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We show how to fully characterize a quantum process in an 
open quantum system. We particularize the procedure to the 
case of a universal two-qubit gate in a quantum computer. 
We illustrate the method with a numerical simulation of a 
quantum gate in the ion trap quantum computer. 
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Recently there has been a growing interest in "quan- 
tum tomography" , i.e. in the complete characterization 
of the state of a quantum system represented by a den- 
sity operator p. Quantum tomography of an unknown 
quantum state (that can be repeatedly prepared) con- 
sists of finding an appropriate sequence of measurements 
which allows one to determine the complete density oper- 
ator p (for experimental implementations and theoretical 
schemes in quantum optics see ]IJ). In this letter we will 
show how to completely characterize a physical process in 
an open quantum system. More specifically, suppose that 
a given quantum dynamics £ transforms input states p- lri 
into output states p ut, i-e. 



Pn 



Pout = £[pin] 



(1) 



with £ a linear mapping. Our aim is to characterize 
the process £, given as a "black box", by a sequence of 
measurements in such a way that it is possible to predict 
what the output state will be for any input state. 

The particular problem that we will analyze after de- 
veloping a general formalism is the characterization of the 
two-bit universal quantum gate for quantum computing 
Pj . A quantum computer consists of n two-level atoms 
with atomic states |0),-, (i — 1, ...,n) representing 
the quantum bits (qubits). States of the quantum com- 
puter are n-atom entangled states in the product Hilbert 
space |V>> € H = n<®W a (i) with « a (») - {|0)i,|l)i}. 
Quantum computations correspond to physical processes 
|Vw) = C^l^in) where a given input state is mapped to 
an output state by a unitary transformation U. This can 
be carried out as a sequence of elementary steps (quan- 
tum gates) involving operations on a few qubits. It has 
been shown that any computation can be decomposed 
into single-bit gates, and a universal two-bit gate which 
involves an entanglement operation on two qubits [Q. 
In reality, due to the presence of decoherence and ex- 
perimental imperfections, these gates (and therefore any 



computation) will not be ideal. In present experiments 
related to quantum computing based on both laser cooled 
trapped ions |3| and atoms in cavities Qj, the difficult 
part is the two-qubit gate, since it requires an interac- 
tion between the two two-level systems via an auxiliary 
system (phonons or photons) which leads to decoherence. 
In view of this fact, we wish to develop a procedure for 
characterizing a two-qubit gate, i.e. characterize a phys- 
ical process £ involving entanglement of two qubits 1 
and 2 in the state space <8> 7^2(2). Below we will 

show how to implement this using only product states as 
inputs, and single qubit measurements on the outputs (as- 
suming that single bit preparations and operations can 
be performed reliably ). We avoid utilizing any interac- 
tion (entanglement) between the qubits which would be 
required to prepare Bell state inputs and perform Bell 
measurements since otherwise the decoherence and er- 
rors induced by the measurement itself would distort the 
characterization of £. Furthermore, we will introduce 
four global parameters to characterize the action of the 
quantum gate: the "Gate Fidelity" (J 7 ), the "Gate Pu- 
rity" (V), the "Quantum Degree" (Q), and the "Entan- 
glement Capability" (£). These four parameters can be 
calculated once the physical process is completely deter- 
mined. To illustrate the procedure of measuring £, we 
will analyze below simulation data from a model of a 
two-bit quantum gate in an ion trap quantum computer 
I- 

To develop the general formalism, let us consider an ex- 
periment in which a quantum system undergoes a phys- 
ical process £. We assume that the system is initially 
prepared in the pure state 
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i=0 



(2) 



where |0), |1), . . . , \N) are orthogonal states spanning the 
Hilbert space of (allowed) input states TL ln with dimen- 
sion N + 1, a subspace of the system Hilbert space 7is- 
We will denote by ps = J2 a u) a \E a )(E a \ the initial state 
of the environment, i.e., of the other degrees of freedom 
that will be coupled to our system. In its most general 
form, the physical process will perform a transformation 
defined by 



1 



\m a ) 



M 

^2\j)\E%), (i = 0,...,N) 

3=0 



(3) 



where in the sum we have taken into account the possible 
presence of other system states that might be populated 
during the interaction (i.e., M > TV in general; see also 
below for the case of a quantum gate). The states \E£j) 
are unnormalized states of the environment. Combining 
(||) and (^|), and tracing over the environment degrees 
of freedom we get the following reduced system density 
operator, 



N 



Pout 



Rr 



(4) 



in the space of output states, H ou t — {|0), . . . , |M}}, and 
where 



M 

RiH= £ b')0"lE w «(^-'l^) (m = o, . . . , n) 

j,j'=0 a 

(5) 

are system operators that do not depend on the ini- 
tial state. From the knowledge of these operators one 
can predict the final density operator p ou t for any input 
state pi n . For a mixed initial state, we diagonalize p m — 

J2 n P n \^i" ){^ia I aim use * ne ^ ac * that * ne ev °lution £ 
is linear, which leads to p out = J2 n Pn £ (\^in^)(Ma D- 
Thus, the problem of fully characterizing the physi- 
cal process £ on the system is reduced to finding the 
(N + l) 2 "transfer operators" R^i in Ti ou t- They fulfill 
Tr{Ej/j} = 5m and (JRf/<)+ = R iV . 

To develop a procedure to measure these operators we 
define two vectorial operators 



Pout : 

R : 



lo {1) 

\ Pout I 
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(6b) 



Here, the components of p ou t are a set of output opera- 
tors corresponding to (N + l) 2 different initial inputs of 

the form (@) with coefficients cf ] [k = 1,...,(N + l) 2 ]. 
All these density operators can be fully characterized us- 
ing standard quantum tomography methods 0. On the 
other hand, the components of the vector R, are defined 
according to R q — Rin, with q = (N + i + In view 
of ((J), these two vectorial operators are related by a set of 
linear equations, written in a matrix form as p ou t = MR, 
where M. is a (c-number) matrix whose elements are de- 
fined as Mk q = c\ k '[c\, ]*. The problem of obtaining 
the transfer operators R thus reduces to finding a set of 
initial states of the system, such that the matrix M. is 
not singular. In that case we will have R = -M -1 p ut, 
which solves the problem. Next, we prove that such an 



invertible matrix M. exists by explicit construction. One 
can simply choose the initial states as follows: 



c[ k) = < 5 ikl if h =k 2 , 

^(^ifci + iSik 2 ) if ki < k 2 , 



(7) 



with k = (N + l)fci + k 2 . 

Let us now apply the above procedure to a universal 
two-qubit gate. Now, the system is composed of two 
two-level subsystems 1 and 2 of levels 10)1,2 and |l)i,2 
each. We can define a set of orthogonal states \i) = 
l K2) 1 (with i = 2i\ + i 2 , and ii,i 2 = 0, 1) and write 
the initial states as in (||) (with dimension N + 1 = 4). 
The quantum tomography of the output states can be 
carried out following the lines proposed by Wootters || : 
one writes the output density operator as 



Pout 



15 
o=0 



(8) 



where A„ 



±q - u qi «y <T 2 2 (q = Aqi + 92), with <r° = 
{l a , a%, &y , o°z\ and a = 1,2 refers to the first and sec- 
ond qubit, respectively. By measuring the observables 
A q , one can determine the coefficients X q , given that 
X q = Tv[p out A q ]/4: 0. Note that all these measurements 
do not require any interaction between the qubits (Bell 
state measurements), that is, for all these measurements 
the two qubits can be measured independently, without 
the application of another two-qubit gate. This is needed 
since otherwise the measurement procedure would in- 
volve errors that could not be separated from the gate 
itself. However, some of the 16 initial states given above 
in order to make up the matrix A4 [cf. Eq. (Q)] are 
entangled states. Their preparation would involve the 
application of a two-qubit gate which would also lead to 
uncontrollable errors. Fortunately, there are other sets of 
initial states that are unentangled for which the matrix 
A4 is invertible. An example are the 16 product states 
|^o)i 1^6)2 (a,b= 1,...,4), where 



|1M = |0>, h/>3> = ^(|0) + |l)) 

m = 11), i^4> = ^(io)+»ii)). 



(9) 



In order to illustrate this procedure, we have studied a 
two-qubit gate in the ion trap quantum computer model 
||. We have considered two ions in a linear ion trap 
interacting with two lasers. Let us denote by \g) n = |0)„ 
and \e) n = \l) n two internal states of the n-th ion, and 
by \e') n an auxiliary internal state. As we have shown in 
Ref . H , the universal two-qubit gate defined by 

ki>i|e 2 >2 -» (-l) eie2 ki)i|e 2 )2, (ei,2 = 0, 1) (10) 

can be implemented in three steps: (i) Apply a it laser 
pulse to the lower motional sideband corresponding to 
the transition \g)x — > |e)i of the first ion; (ii) Apply a 2ir 
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laser pulse to the lower motional sideband of the transi- 
tion | (7)2 — > |e'/2 of the second ion; (hi) as (i). By lower 
motional sideband we mean that the laser frequency has 
to be equal to the corresponding internal transition fre- 
quency minus the trap frequency, in order to excite a 
center of mass phonon only. The interaction of the two 
ions and the laser is given by the following Hamiltonian 
I 

H = -Ai|e)n(e| - A 2 |e') 2 2<e'| + va\ m a cm + \f3va\a x 
+ -^[|e)n(ff|e" 4 ' , ™ (a ™ +a ™ ) e- 4 ''' (Q ' +a ' t) + H.c] 

Here, Ai^ and 0^2 are the laser detunings and Rabi 
frequencies of the laser acting on each ion, respectively. 
The operators a and a* are annihilation and creation op- 
erators of the center of mass (cm) and relative (r) mo- 
tion mode, 77 is the corresponding Lamb-Dicke parame- 
ter, and v is the trap frequency. 

We have calculated numerically the evolution given by 
this Hamiltonian. Following the steps mentioned above, 
and for Rabi frequencies much smaller than the trap fre- 
quency this evolution corresponds to the gate jlO| ) ||. 
For finite Rabi frequencies, however, the result will not 
be ideal. After the gate operation, there will remain some 
population in the phonon modes, which would lead to 
decoherence. Moreover, there may remain some popu- 
lation in the auxiliary state |e')2, and therefore in this 
case we have M = 4 > N = 3. With the numeri- 
cal calculation we have simulated the measurement of 
the operators Ri>i, that fully characterize the evolution 
process. We emphasize that all the information regard- 
ing this (non-ideal) gate is contained in these operators. 
In Fig. 1 we have compared the ideal case of a perfect 
gate [Fig. 1(a)] with the simulation results with realis- 
tic parameters [Figs, l(b-d)]. We have plotted the ma- 
trix elements of 16 operators R^i sorted according to 
£n, m = with n — 4i + j and m = 4i' + j' . We 

have chosen a set of parameters close to those planned in 
experiments: rj cm = r) r = 0.5, Ai = A2 = —v, and Rabi 
frequencies Vt\ — VI2 — 0.1,0.2z/ and 0.5 [Figs.l(b,c,d)] 
|1. As it is shown, for moderately small Rabi frequen- 
cies the simulated results almost coincide with the ideal 
ones. 

In order to elucidate to what extent a two-qubit gate 
implemented experimentally in a particular model of 
quantum computer approaches to the ideal one, we define 
a parameter, the "Gate Fidelity", as 

T= (^inl^Poutt/^in), 

where the overline indicates average over all possible in- 
put states l^in), and U is the unitary operator corre- 
sponding to the ideal gate. This parameter can be cal- 
culated once the full characterization of the gate is per- 
formed using 



i=0 vjtj 

where Fj,j = Ri<iU\j) . Obviously, a gate fidelity 

close to one indicates that the gate was carried out almost 
ideally. 

In a similar way we can define the "Gate Purity" 
V = Tr{(p out ) 2 } that reflects the effects of decoherence 
on the gate. V close to one indicates that the effects of 
decoherence are negligible. It can be shown that 

1 3 1 
V = - Tr{(i?») 2 } + ^4 ^{RiiRi'V + Ri'iRii'}- 

i=0 i^j 

In addition, the "Quantum Degree of the Gate" Q is 
defined as the maximum value of the overlap between 
all possible output states that are obtained starting from 
an unentangled state and all the maximally entangled 
states, i.e. 

Q = max (* mc |/0out|^me), 

Pout,|*mc> 

where p on t denote the output states corresponding to un- 
entangled input states l^m) = |V'a)i|^'i>)2, and |^ roe ) is 
a maximum entangled state As it has been shown, 
when the overlap between a density operator and a max- 
imally entangled state is larger than (2 + 3v / 2)/8 ~ 0.78, 
Clauser-Horne-Shimony-Holt (CHSH) inequalities are 
violated fPlJ| , Finally, another useful parameter is the 
"Entanglement Capability" £ pQ ], given as the smallest 
eigenvalue of the partial transposed density matrix /5 ou t , 
for unentangled inputs states. As it has been recently 
shown |]lo|| , the negativity of this quantity is a necessary 
and sufficient condition for non-separability of density 
operators of two spin-1/2 systems. These quantities can 
be calculated numerically starting from the gate opera- 
tors Ri'i with a maximization/minimization procedure. 

In Fig. 2 we have plotted these four parameters as a 
function of the Rabi frequency = Ox — ^2 for three dif- 
ferent Lamb-Dicke parameters and the same parameters 
as in Fig. 1. As expected, the best results are obtained 
for small Rabi frequencies. This is due to the fact that 
transitions to undesired levels are suppressed. On the 
other hand, a Lamb-Dicke parameter close to one also 
improves the results. This is due to the fact that in con- 
trast to Ref. H we are considering here traveling wave 
excitation, in which transitions that do not change the 
phonon number can be excited by the laser. The ratio 
between the effective Rabi frequencies of the lower side- 
band and these unwanted ones is proportional to rj in the 
Lamb-Dicke limit. The figures also indicate that the T , 
V, <2, and £ decay in a different manner when the pa- 
rameters deviate from the conditions of operation of the 
ion trap quantum computer. 

Finally, we wish to mention another application of this 
procedure of characterization of a quantum physical pro- 
cess. Let us consider a quantum system coupled to a 
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Markovian reservoir. In the case that the Hamiltonian is 
time independent, one can describe the evolution of the 
reduced system density operator p in terms of a master 
equation of the form p = Cp, where C is a Liouvillian su- 
peroperator. The formal solution to this equation after 
a time t is p ou t = e ct pi n . The full characterization of the 
physical process in this case would allow one to "mea- 
sure" the Liouvillian C, and therefore to determine the 
master equation fulfilled by the system. Moreover, choos- 
ing different interaction times one could check whether a 
given process is Markovian or not. 

In conclusion, we have shown how to perform the full 
characterization of a physical process. This requires the 
preparation of various initial states, and the quantum to- 
mography of the corresponding output states. We have 
illustrated these results for the case of a quantum gate, 
and have defined four parameters that give the quality 
of the gate. In addition, we have presented numerical 
simulations for the ion trap quantum computer. The 
complete characterization of the quantum gates is rele- 
vant both from an experimental point of view to evaluate 
implementations of quantum gates in the laboratory, as 
well as a theoretical tool to compare the expected per- 
formance of specific quantum computer model systems. 
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(b,c,d) 
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FIG. 2. Fidelity (a), Purity (b), Quantum Degree (c) and 
Entanglement Capability (d) for a two-qubit gate in the ion trap 
quantum computer as a function of the Rabi frequency. The 
values of the Lamb-Dicke parameter are indicated. All the other 
parameters are as in Fig. 1 
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